ABSTRACT 

ireless  mesh  networks  are  sys¬ 
tems  of  interconnected  wireless 
access  points  that  provide  digital 
services  to  client  devices  via  radio  trans¬ 
mission.  We  consider  the  challenges  of  a 
communications  planner  who  must  quickly 
design  a  wireless  mesh  network,  as  might 
be  expected  during  combat  operations  or 
in  support  of  humanitarian  assistance  and 
disaster  relief  operations.  We  seek  a  network 
that  maximizes  client  coverage  area  subject 
to  constraints  on  network  service,  the  tech¬ 
nical  characteristics  of  the  available  access 
points,  and  radio  propagation  over  terrain. 
We  create  a  nondifferentiable,  nonconvex, 
nonlinear  optimization  problem  and  use  a 
sampling  algorithm  to  quickly  find  good 
solutions.  We  validate  our  formulation  and 
solutions  via  numerical  experiments  and 
several  field  tests,  and  we  demonstrate  that 
our  technique  can  generate  network  topol¬ 
ogies  capable  of  functioning  in  real-world 
scenarios. 

We  construct  a  corresponding  decision 
support  tool  that  allows  a  communications 
planner  to  design  working  wireless  mesh 
network  topologies  quickly,  with  no  guess¬ 
work,  and  requiring  very  little  expertise. 
The  tool  runs  on  a  laptop,  supports  virtu¬ 
ally  any  type  of  access  point,  uses  terrain 
information  freely  downloadable  from  the 
Internet,  and  does  not  require  any  addi¬ 
tional  software  or  solver  licenses. 

INTRODUCTION 
Description  of  Problem 

A  wireless  mesh  network  (WMN)  is 
a  communications  network  of  fixed  access 
points  (APs)  that  exchange  electronic  mes¬ 
sages  via  radio  transmission  to  and  from  cli¬ 
ent  devices  (such  as  computers,  sensors,  or 
mobile  devices).  The  fixed  position  of  its 
APs  differentiates  a  WMN  from  a  so-called 
ad  hoc  network,  where  the  APs  can  be  con¬ 
stantly  moving  (Zhang  et  al.  2006,  p.  565). 

Military  and  civilian  organizations  can 
benefit  from  the  inherent  advantages  of 
WMNs.  During  combat  operations,  WMNs 
can  quickly  and  securely  relay  time-critical 
information  such  as  intelligence  reports, 
tactical  orders,  and  sensor  readings  to  sep¬ 
arated  small  units.  During  humanitarian 


assistance  and  disaster  relief  (HA /DR)  op¬ 
erations,  WMNs  can  provide  maps,  floor- 
plans,  video  surveillance,  emergency  aid 
requests,  and  other  critical  information  to 
first  responders. 

In  the  system  under  consideration,  each 
AP  uses  two  separate  and  configurable  ra¬ 
dio  devices  to  create  a  WMN  with  two 
levels  of  connectivity.  The  first  level  of  con¬ 
nectivity  supports  client-to-AP  communi¬ 
cation  within  a  client  coverage  area,  whereas 
the  second  level  consists  of  a  backhaul  radio 
network  that  routes  traffic  between  APs.  Cli¬ 
ent  devices  may  roam  within  the  coverage 
area  and  communicate  with  one  another 
or  to  an  outside  network,  such  as  the  Inter¬ 
net,  through  a  gateway  (e.g.,  via  a  satellite 
uplink).  We  assume  that  wireless  APs  alone 
provide  traffic  routing  services,  thus  form¬ 
ing  an  infrastructure  mesh  type  of  WMN 
(see  Nicholas  2009,  pp.  2-4).  Additionally, 
we  assume  client  devices  communicate  with 
a  single  AP  at  a  time,  and  hence  do  not  serve 
as  intermediate  relay  points  (Zhang  et  al. 
2006,  pp.  564-567). 

The  physical  topology  of  a  WMN,  as  de¬ 
fined  by  the  locations  of  the  wireless  APs,  is 
critical  to  its  performance.  We  consider  the 
challenge  of  a  communications  officer  or 
a  network  designer  who  must  choose  the  loca¬ 
tions  and  configurations  of  wireless  APs  to 
provide  service  to  clients  in  desired  cover¬ 
age  areas  while  meeting  restrictions  on  the 
quantity,  placement,  and  characteristics  of 
the  APs  and  also  satisfying  any  require¬ 
ments  for  coverage,  bandwidth,  and  other 
service  standards.  The  designer  must  also 
consider  the  effects  of  terrain  and  other  as¬ 
pects  of  the  operating  environment  on  ra¬ 
dio  wave  propagation.  Because  combat 
and  HA /DR  operations  are  time-sensitive, 
the  designer  must  build  the  WMN  quickly 
and  with  as  little  guesswork  as  possible. 

We  focus  on  designing  a  WMN  topol¬ 
ogy  that  maximizes  client  coverage  while 
considering  network  flow  requirements. 
There  is  a  fundamental  tension  when  de¬ 
signing  such  a  WMN:  maximizing  client 
coverage  tends  to  place  APs  far  apart, 
whereas  the  requirements  for  backhaul  net¬ 
work  throughput  tend  to  keep  APs  relatively 
close.  The  challenge  is  to  balance  this  trade¬ 
off  for  a  fixed  number  of  APs  while  re¬ 
specting  their  technological  details  and  the 
physics  of  radio  propagation  over  terrain. 

We  evaluate  the  performance  of  a  given 
WMN  topology  in  two  steps.  First,  we 
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calculate  the  value  of  client  coverage  as  dictated 
by  specific  AP  locations  and  their  configura¬ 
tion,  ground  terrain,  and  environmental  data. 
We  then  calculate  the  value  of  network  flow 
according  to  the  Simultaneous  Routing  and 
Resource  Allocation  (SRRA)  techniques  of  Xiao 
et  al.  (2004).  We  combine  these  two  subprob¬ 
lems  to  formulate  SRRA+C,  a  nondifferentiable, 
nonconvex,  nonlinear  optimization  problem.  We 
implement  the  Dividing  RECTangles  (DIRECT) 
sampling  algorithm  of  Jones  et  al.  (1993)  to 
quickly  find  good  solutions  to  SRRA+C. 

To  validate  the  quality  of  the  network  topol¬ 
ogies  generated  by  our  technique,  we  conduct 
numerical  experiments  and  field  tests  using  com¬ 
mercial  equipment.  As  an  aid  to  the  network 
designer,  we  create  a  customized  graphical  de¬ 
cision  support  tool  to  solve  the  SRRA+C  prob¬ 
lem.  The  standalone  tool  reads  and  graphically 
displays  digital  terrain  elevation  information, 
obtains  its  best  solution  to  the  SRRA+C  prob¬ 
lem,  and  displays  the  resulting  network  and 
client  coverage  regions.  We  demonstrate  that 
our  technique  can  quickly  create  WMN  topolo¬ 
gies  that  function  in  realistic  scenarios,  and  that 
our  decision  support  tool  can  assist  communi¬ 
cations  officers  or  network  designers  in  build¬ 
ing  WMNs  in  support  of  combat  or  HA/DR 
operations. 

Previous  Work 

He  et  al.  (2004)  use  the  DIRECT  sampling  al¬ 
gorithm  of  Jones  et  al.  (1993)  to  find  good  AP 
placements  in  indoor  wireless  networks,  using 
very  accurate  (and  computationally  expensive) 
ray-tracing  techniques  that  predict  radio  prop¬ 
agation.  We  build  on  this  general  idea  to  cal¬ 
culate  the  value  of  client  coverage,  and  use 
the  DIRECT  algorithm  similarly  for  outdoor 
environments. 

Xiao  et  al.  (2004)  solve  the  SRRA  problem 
via  dual  decomposition  to  identify  optimal  traf¬ 
fic  routes  and  allocation  of  AP  transmission 
power  within  a  wireless  network.  As  noted 
above,  we  use  SRRA  to  calculate  the  value  of 
network  flow. 

Shankar  (2008)  uses  the  SRRA  formulation 
to  determine  network  flow  among  preposi¬ 
tioned  wireless  nodes,  and  then  adopts  the 
attacker-defender  techniques  of  Brown  et  al. 


(2006)  to  calculate  optimal  jammer  locations 
that  maximally  disrupt  network  flow.  Our 
work  complements  Shankar's  research,  pro¬ 
viding  the  initial  design  of  a  WMN. 

To  our  knowledge,  ours  is  the  first  tech¬ 
nique  for  designing  WMN  topologies  that 
maximize  client  coverage  while  considering 
backhaul  network  service  requirements  on  real 
terrain.  See  Nicholas  (2009,  pp.  6-9)  for  a  de¬ 
tailed  literature  review. 

This  paper  is  organized  as  follows.  In  the 
next  section,  we  describe  each  element  of  the 
SRRA+C  formulation  along  with  our  technique 
to  solve  it.  We  then  consider  several  notional 
and  real  network  design  problems,  contrast  our 
solution  to  what  we  obtain  from  brute  force 
enumeration,  and  briefly  summarize  the  results 
of  our  field  tests.  We  conclude  with  suggestions 
for  follow-on  research. 


MODEL  FORMULATION 

Our  goal  is  to  position  APs  in  locations  that 
maximize  client  coverage,  subject  to  restric¬ 
tions  on  network  service,  AP  characteristics 
and  placement,  and  radio  propagation  over 
terrain.  As  noted  earlier,  there  is  an  inherent 
tension  between  maximizing  client  coverage 
and  network  traffic  flow.  The  crux  of  the 
SRRA+C  formulation  is  to  capture  and  quan¬ 
tify  this  tension. 

We  represent  each  AP  as  a  node  in  a  net¬ 
work.  Let  N  denote  the  set  of  all  AP  nodes, 
indexed  by  i  =  1,  2,  ...,  n,  where  n  =  \N\.  Let 
A  =  (A, .  A2  ?=  A„)  represent  the  locations  of  the 
nodes.  Although  the  mathematics  support 
any  coordinate  system  to  represent  location, 
in  our  implementation  each  node  location  A, 
is  itself  a  two-dimensional  coordinate  repre¬ 
senting  the  northing  and  easting  for  node  i. 
We  assume  that  node  locations  are  fixed  in  the 
sense  that  nodes  will  remain  in  position  once 
placed.  In  what  follows,  we  use  the  terms  AP 
and  node  interchangeably. 

We  define  the  operating  region  as  the  topo¬ 
graphic  area  where  an  AP  may  be  physically 
located.  We  partition  the  operating  region 
into  a  set  of  discrete  regions  R,  indexed  by 
r  =  1,  2,  ...  \R\.  Although  our  formulation  al¬ 
lows  the  use  of  any  discretization  scheme. 
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our  implementation  assumes  rectangular  re¬ 
gions  arranged  in  a  grid  (see  Figure  1). 

Each  AP  has  two  radio  devices.  The  first 
provides  client  coverage  in  the  immediate  vi¬ 
cinity  of  the  AP;  the  specific  coverage  obtained 
(e.g.,  shaded  area  in  Figure  1)  depends  on  sev¬ 
eral  factors  including  the  local  terrain.  The 
second  radio  creates  a  backhaul  network  be¬ 
tween  nodes  (e.g.,  dashed  lines  in  Figure  1). 
Let  AEN  XN  denote  the  set  of  directed  back¬ 
haul  arcs  between  nodes,  with  flow  along  arc 
(i,j)  e  A  representing  the  directed  transmission 
from  node  i  to  node  j.  In  general,  the  wireless 
nature  of  this  network  means  that  transmis¬ 
sion  along  any  of  the  backhaul  arcs  is  possible, 
although  distance,  terrain,  and  background  in¬ 
terference  will  dictate  which  of  these  arcs  has 
the  transmission  capacity  to  be  useful.  Follow¬ 
ing  Xiao  et  al.  (2004),  we  assume  APs  are  not 
subject  to  self-jamming  or  interference  from 
other  APs. 

We  identify  a  single  node  d  £  N  as  the  head¬ 
quarters  (HQ)  node  (see  Figure  1),  and  assume 


the  network  designer  predefines  the  location 
of  this  node.  (In  what  follows,  we  will  assume 
that  node  d  is  always  the  first  node  in  the  set  N.) 
We  assume  the  vast  majority  of  client  network 
traffic  will  be  directed  to  or  received  from  the 
HQ  node,  as  this  location  will  connect  to  the 
Internet  or  other  outside  network,  as  well  as 
host  email,  domain,  and  storage  servers.  Hence, 
we  optimize  our  network  for  traffic  flow  from 
client  service  areas  to  the  HQ  node,  though 
our  formulation  is  more  general  and  allows 
any  number  of  APs  to  serve  as  the  destination 
for  network  traffic. 

Calculating  Client  Coverage 

The  client  coverage  provided  by  a  particu¬ 
lar  WMN  topology  is  a  function  of  its  AP  loca¬ 
tions.  We  calculate  the  received  signal  strength 
(RSS)  pir  (measured  in  dBm)  from  the  trans¬ 
mitter  at  node  i  to  the  center  of  region  r  using 
the  standard  link  budget  formula  (Olexa  2005, 
p.  79): 


A 

r 

✓ 

J 

Q 

1 

✓ 

t 

\ 

\ 

X 

\ 

X 

4 

/ 

✓ 

\ 

\ 

V 

' 

* 

X 

✓ 

\ 

X 

r 

hi 

\ 

\ 

✓ 

f 

/s 

X 

\ 

\ 

/ 

/ 

\ 

X 

> 

\ 

J 

L' 

\ 

X 

> 

Q 

*"*  «» 

—  

X 

X 

-- 

X 

-y 

] 

□ 

Figure  1.  Discretized  operating  region  and  wireless  mesh  network  (WMN).  Circles  denote  the  location  of 
access  points  (APs),  shaded  regions  denote  the  areas  with  sufficient  client  coverage  (i.e.,  zero  coverage 
shortfall),  and  dashed  lines  denote  the  backhaul  network.  In  this  example,  there  is  a  single  headquarters 
(HQ)  node  (labeled  d)  and  four  APs  (labeled  2,  3,  4,  5).  The  locations  of  the  APs  are  represented  by  the  vector 
A  =  (Arf,  A2,  A3,  A4,A5). 
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pir  =  power j  +  gaitii  -  lossi  -  losspath 

—  lossmisc  +  gainr  —  los  sr  (1) 

where  power,  is  transmitted  power  (in  dBm) 
from  node  i;  gain;  and  gain,,  are  respectively 
the  gains  (in  dBi)  of  the  transmitter  at  node  i 
and  receiver  in  region  r;  loss,  and  lossr  are  re¬ 
spectively  the  losses  (in  dB)  from  cables,  con¬ 
nectors,  etc.,  of  the  transmitter  and  receiver; 
losspath  is  total  path  loss  (in  dB)  from  A,  to  the 
center  of  region  r;  and  /m,-sc  is  miscellaneous  loss 
(in  dB),  such  as  fade  margin.  All  of  the  terms 
in  Equation  (1)  are  input  data,  determined  by 
the  equipment  and  environment,  except  for  the 
total  path  loss  losspath,  which  depends  on  the 
position  of  the  transmitter  and  receiver  in 
question. 

Our  formulation  allows  any  model  for  com¬ 
puting  losSpath,  including  the  Irregular  Terrain 
Model  (ITM)  (Longley  and  Rice  1968)  and 
Hata-COST  231  (COST  1999).  Our  preferred 
model  is  the  Terrain  Integrated  Rough  Earth 
Model  (TIREM)  of  Alion  Science  &  Technology 
Corporation  (http:/ /www.alionscience.com).  This 
model  samples  terrain  elevation  to  compute  path 
loss,  and  considers  the  effects  of  free  space  loss, 
diffraction  around  obstacles,  atmospheric  ab¬ 
sorption  and  reflection,  and  other  factors. 

We  adopt  and  modify  He  et  al.'s  (2004)  con¬ 
cept  of  power  coverage  to  quantify  the  value  of 
client  coverage.  We  specify  for  each  region 
r  £  R  a  minimum  coverage  threshold  t,  in  dBm; 
any  received  signal  p,r  above  this  threshold 
from  AP  i  to  region  r  qualifies  as  adequate  cli¬ 
ent  coverage.  The  difference  between  p,>  and 
Tr  represents  a  quantity  we  define  as  coverage 
shortfall  at  region  r  from  node  i.  Summarizing, 
for  each  location  r  £  R  and  each  node  i  £  N, 
we  have 

(Coverage  Shortfall),-,  =  (rr  —  pir)  + 

where  ()+  denotes  the  projection  onto  the  non¬ 
negative  real  line.  Because  a  positive  differ¬ 
ence  represents  inadequate  client  coverage, 
we  wish  to  minimize  this  quantity.  We  need 
consider  only  the  minimum  coverage  shortfall 
from  each  node  i,  as  we  assume  each  client 
device  can  connect  to  only  one  AP.  We  sum 
over  all  r  £  R  to  calculate  total  coverage  short¬ 
fall,  denoted  Z coverage * 


Z coverage  )  =  ( Total  Coverage  Shortfall) 


The  total  coverage  shortfall  is  a  function  of 
node  locations  A.  By  allowing  only  positive 
terms,  we  disallow  the  benefit  of  transmitting 
excessive  power  to  any  given  coverage  location. 
This  method  also  limits  the  amount  any  cover¬ 
age  location  can  be  penalized:  no  more  than 
Tr.  Note  it  is  possible  to  increase  the  relative  im¬ 
portance  of  any  region  r  by  multiplying  the  re¬ 
gion's  coverage  shortfall  by  a  positive  scalar. 

Calculating  Network  Flow 

Given  fixed  AP  locations  A,,  i  £  N,  we  use 
the  Shannon  capacity  formula  (1949)  to  calcu¬ 
late  the  transmission  capacity  along  each  arc 
(i.j)  £  A.  This  formula  establishes  a  theoretical 
upper  bound  on  transmission  capacity  in  bits 
per  second  (bps).  Following  Xiao  et  al.  (2004), 
the  capacity  from  node  i  to  node  j  is: 

(Capacity),,  -  bandwidth  log2(l  +  ,j) 

mi)  ft  A  (3) 

where  bandwidth  is  channel  bandwidth  in  Hertz; 
gain,)  is  the  sum  of  the  antilog  gain  terms  (gainj 
and  gainj  from  Equation  1);  noise p  is  the  back¬ 
ground  noise  power  in  watts  from  node  i  to  node 
j;  and  loss,,  is  the  sum  of  the  antilog  loss  terms 
(loss,-,  losSj,  losSpati„  and  Zossm,sc  from  Equation  1). 
Again,  these  input  data  depend  on  the  place¬ 
ment  of  APs,  which  we  assume  are  at  known, 
fixed  locations.  In  theory,  we  assume  that  all 
backhaul  arcs  are  possible  (i.e..  A  —  N  X  N). 

We  assume  each  AP  has  limited  total  trans¬ 
mission  power  denoted  p,  (in  watts),  and  we 
define  P,;  to  be  the  amount  of  p,  used  to  trans¬ 
mit  from  i  to  j.  Thus,  each  AP  is  additionally 
constrained  by 

E  pn^Vi-  (4) 

Here,  P,;  is  a  decision  variable  representing  the 
AP-to-AP  transmission  power  from  node  i  to 
node  j,  whereas  the  transmission  power  for 
AP-to-client  communication  power ,  is  a  (con¬ 
stant)  input  parameter. 
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We  measure  each  individual  traffic  flow 
in  bps.  We  adopt  the  approach  of  Xiao  et  al. 
(2004)  to  quantify  the  value  of  total  network 
flow  according  to  a  log-utility  function  that 
places  a  zero  value  on  unit  flow,  positive  values 
on  flows  greater  than  one,  and  negative  values 
on  flows  less  than  one.  We  note  that  a  zero 
flow  has  an  infinite  penalty  and  therefore  there 
is  strong  incentive  to  ensure  that  each  origin- 
destination  pair  receives  some  flow.  Defining 
Sf  to  be  the  total  flow  originating  at  node  i  and 
destined  for  node  d,  we  have 

(Utility  of  Total  Network  Flow)  log2(Sf) . 

d  i^d 

(5) 

Collectively,  we  obtain  our  version  of  the  Xiao 
et  al.  (2004)  SRRA  problem  for  given  AP  loca¬ 
tions  A  as  shown  in  Figure  2. 


Given  AP  locations  A,  this  is  a  multicom¬ 
modity  network  flow  problem.  The  objective 
function  (SO)  maximizes  the  total  utility  of 
traffic  flow  between  each  origin-destination 
pair.  Constraints  (SI)  ensure  balance  of  flow 
at  each  node.  Constraints  (S2)  define  the  total 
flow  along  any  arc  as  the  sum  of  all  traffic 
flows  along  that  arc.  Constraints  (S3)  ensure 
that  total  flow  along  any  arc  is  less  than 
or  equal  to  its  capacity.  Constraints  (S4)  en¬ 
force  a  budget  on  the  total  transmission 
power  at  each  AP.  Constraints  (S5)-(S8)  ensure 
nonnegativity. 

Xiao  et  al.  (2004)  observe  that  the  SRRA 
problem  can  be  solved  via  dual  decomposi¬ 
tion  because  of  its  layered  structure.  Specifi¬ 
cally,  by  introducing  the  Lagrange  multipliers 
ayV(zj)  e  A  for  the  capacity  constraint  (S3), 
we  obtain  the  partial  Lagrangian 


Index  Use 

ieN 

(UJ)*A 

deD^N 


node  (alias  j) 
directed  arc  (link) 
destination  nodes 


Input  Data 

X, 


Pi 

bandwidth 


the  locations  of  nodes.  A,  =  {A,-,  ieN} 
maximum  total  transmission  power  per  node,  ieN 
channel  bandwidth 


[none] 


[watts] 

[hertz] 


Calculated  Data 

gairij 

lossy 

noisej 


sum  of  antilog  gain  terms  from  ieN  to  j  e  N 
sum  of  antilog  loss  terms  from  ieN  to  j  e  N 
background  noise  power  from  ieN  to  j  eN 


[none] 

[none] 

[watts] 


Decision  Variables 


sf 

total  flow  of  traffic  from  origin 

ieN  to  destination  d eD 

[bps] 

N 

traffic  flow  along  arc  (i,  j )  e  A  to  destination  d  eD 

[bps] 

Ta 

total  flow  along  arc  (i,  j)  e  A 

[bps] 

pu 

total  transmission  power  along  < 

arc  ( i,j)eA 

[watts] 

Formulation 

W  =  may  E  X  lo&2  (si  ) 

*’  ’  ’  d  i*d 

(SO) 

s.t.  X  F*- 

XJ,i)eA 

y  Fd=sd 

La  y  J 

i:(i,j)eA 

V /  e  N,  Vrf  e  D 

(SI) 

V 

W' 

II 

V(i,j)eA 

(S2) 

Tu  -  bandwidth  log  J  1  +  — — — - —  R.  <0 
"  1  noisejlosSij  J 

V(i,j)eA 

(S3) 

I  r,s 

P, 

Vi  g  N 

(S4) 

Si  >0 

i^d 

(S5) 

f;>  o 

V(i,j)eA,  Vd  eD 

(SO) 

Ttj>0 

V(i,j)eA 

(S7) 

P„  ao 

V(i,j)eA 

(S8) 

Figure  2.  Formulation  SRRA. 
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L(S,F,T,P,a)  =  J211  l°g2(S/) 

d  i^d 

—  ad  Tij  —  bandwidth 

M)  eA  ' 


log2 


8ainij 

noiseijlossij 


(6) 


The  objective  function  of  the  dual  problem  is: 


V(a)  =  max  L(S,F,T,P,  a).  (7) 

One  immediate  observation  is  that  this  dual 
function  can  be  evaluated  separately  in  the  net¬ 
work  flow  variables  S,  F,  and  T  and  the  com¬ 
munications  variable  P.  Thus,  the  problem  of 
evaluating  the  dual  function  separates  into 
a  network  flow  (net)  subproblem  and  a  resource  allo¬ 
cation  (comm)  subproblem,  that  is: 


V(a)  =  Vnet(a)  +  Vcomm(a)  (8) 

where 


)“^EElos 2(s?)- 

aijTij 

d 

i,i^d 

(v)6A 

s.t. 

E 

E  d 

=  Sd 

Vj  €  N,  \/d  €  D 

i: 

VI)  6+  i 

P,j)eA 

II 

-M 

^Tl 

d 

V 

V(i;)eA 

S?>  0 

i  ^  d 

O 

Al 

"a 

Us 

V(i,;)  eA, 

VdeD 

0 

V(i,;)eA 

V commipt) 

av 

bandwidth 

(v)eA 

l0g2 

(i+  *7  p,) 

\  noise ij  lossij  ‘ ) 

s.t. 

E  J> 

ij^Vi 

VieN 

J:(V)6A 

Pij^  0 

V(i,/)  < 

EA. 

The  Lagrange  dual  problem  is  thus: 


min  V(a)  =  Vnet(a)  +  Vcomm(a)  (9) 

a 

s.t.  a  >  0. 


The  dual  function  17(a)  is  always  convex.  Xiao 
et  al.  (2004)  assume  that  a  feasible  solution 
(S,F,T,P)  exists  such  that  the  nonlinear  capacity 
constraints  hold  with  strict  inequality  (known 
as  Slater's  condition,  see  Boyd  and  Vandenberghe 
2004,  Sec.  5.2),  and  therefore  conclude  strong 
duality  holds  and  the  optimal  value  of  this  dual 
problem  is  equal  to  the  optimal  value  of  the  pri¬ 
mal  problem.  However,  they  also  note  that  the  ob¬ 
jective  function  of  the  primal  problem  is  not 
strictly  concave  in  the  variables  F  and  T,  and  thus 
the  dual  function  is  only  piecewise  differentiable. 
As  a  result,  the  dual  problem  is  a  nondifferentiable 
convex  optimization  problem,  to  which  they  ap¬ 
ply  the  subgradient  method  to  obtain  a  solution. 

We  also  apply  the  subgradient  method  to 
solve  this  problem,  running  it  only  a  fixed 
number  of  iterations  (typically  500)  to  obtain  an 
approximate  solution.  Each  iteration  of  the  sub¬ 
gradient  method  might  not  necessarily  improve 
the  dual  objective  value,  but  each  iteration  re¬ 
duces  the  optimality  gap  (Bertsekas  1999,  p.  621). 

As  noted  by  Xiao  et  al.  (2004),  both  subprob¬ 
lems  are  convex  optimization  problems  with 
special  structure  lending  themselves  to  very  ef¬ 
ficient  computational  techniques.  We  solve  the 
network  flow  subproblem  as  a  multicommodity 
network  flow  problem,  and  use  an  algorithm 
described  by  Luss  and  Gupta  (1975)  to  solve 
the  resource  allocation  subproblem.  Further  de¬ 
tails  of  the  subgradient  algorithm  and  its  imple¬ 
mentation  to  solve  this  problem  are  available  in 
Nicholas  (2009,  pp.  35-45).  We  use  the  approxi¬ 
mate  solution  to  the  SRRA  problem  as  a  quanti¬ 
fication  of  the  value  of  network  flow. 

SRRA+C  Formulation 

The  overall  SRRA+C  problem  minimizes 
Z(A)by  choice  of  AP  locations  A„  Vi  e  N.  We  cal¬ 
culate  the  overall  performance  of  a  WMN  as: 

Z(A)  =  Zcoverage  (A )  -  wZfhw (A ) .  (10) 

That  is,  the  overall  objective  function  is  a  linear 
combination  of  client  coverage  (calculated  as 
client  coverage  shortfall)  and  network  flow. 
The  objective  function  is  in  units  of  dBM,  al¬ 
though  the  objective  value  has  no  direct  practi¬ 
cal  interpretation.  Rather,  the  objective  value 
serves  as  a  relative  method  of  comparing  differ¬ 
ent  network  topologies. 
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We  use  if  as  a  positive  scalar  (in  units  of 
dBm/log2  bps)  representing  the  relative  weight 
placed  on  network  flow.  Larger  values  of  w  in¬ 
crease  the  value  of  network  flow,  and  in  general 
increase  the  appeal  of  more  compact  network 
topologies.  We  typically  use  w  =  1,  meaning  that 
we  weight  these  terms  equally.  See  Nicholas 
(2009,  pp. 78-80)  for  a  detailed  sensitivity  analy¬ 
sis  of  w. 

Solving  the  SRRA+C  Problem 

One  method  of  solving  the  SRRA+C  prob¬ 
lem  is  simply  to  restrict  APs  to  a  discrete  set  of 
fixed  locations  (e.g.,  the  center  of  each  location 
r  €  R)  and  then  try  all  possible  network  solu¬ 
tions  via  complete  enumeration.  This  requires 
that  we  calculate  the  overall  objective  value  for 
each  unique  topology  and  keep  track  of  the 
best  one(s).  Recall  that  there  are  |R|  operating  re¬ 
gions  and  n  APs,  with  the  position  of  the  first  AP 
(i.e.,  the  HQ  node)  already  fixed.  Hence,  the  to¬ 
tal  number  of  unique  topologies  is 

and  the  solution  space  grows  exponentially  in 
both  |R|  and  n.  Clearly,  we  need  a  more  efficient 


method  of  solving  SRRA+C  if  it  is  to  be  used  for 
fast  design  of  WMNs.  We  show  how  to  use  the 
Dividing  RECTangles  (DIRECT)  algorithm  of 
Jones  et  al.  (1993)  to  meet  this  need. 

Let  A1'  =  (AjjAj,  A*)  denote  the  positions 
of  the  APs  in  iteration  k  of  the  algorithm.  We 
follow  the  algorithm  in  Figure  3  to  iteratively 
select  A/c  for  increasing  k  until  we  reach  a  speci¬ 
fied  stopping  criterion  (number  of  iterations). 

This  iterative  scheme  is  a  modified  version 
of  DIRECT,  which  is  a  sampling  optimization 
algorithm  based  on  Lipschitzian  optimization 
(Horst  and  Hoang  1996,  pp.  43-46)  that  iteratively 
divides  the  multidimensional  solution  space 
of  AP  locations  into  smaller  hyper-rectangles 
(hence  the  name  Dividing  RECTangles).  Un¬ 
like  Lipschitzian  optimization,  DIRECT  does  not 
require  a  priori  specification  of  the  Lipschitz  con¬ 
stant,  nor  knowledge  of  the  objective  function 
gradient,  which  makes  it  appealing  in  solving 
this  non-differentiable,  non-convex,  nonlinear 
problem. 

In  this  implementation,  DIRECT  samples 
candidate  solutions  within  a  unit  hypercube 
defined  by  the  multidimensional  SRRA+C  so¬ 
lution  space.  The  dimensionality  of  this  space 


Input:  Full  SRRA  problem  data  (number  and  operating  characteristics  of  APs,  HQ  node  location, 
elevation  and  coverage  requirements  for  each  r  e  R)  and  desired  number  of  iterations 
max -iterations. 

Output:  Best  estimate  of  optimal  AP  locations  A*  =  A*  J  and  overall  objective  value,  Z^A* ) . 

begin 

Store  map  data 

Initialize  starting  AP  locations  A1  at  center  of  unit  hypercube;  calculate  Z(Al) 

Initialize  k  <—  2 

while  (k  <  max -iterations)  do 

Select  candidate  AP  locations  Ak  based  on  previous  objective  values  Z(A1),^±,Z(Ak~I) 
and  hyper-rectangles  for  previous  solutions  A\...,Ak~l 
Calculate  coverage  shortfall  value  Zcoverage  (  Ak )  according  to  Equation  (2) 

Solve  formulation  SRRA  and  calculate  network  flow  values  Zflow  ( Ak ) 

Calculate  overall  objective  value  Z  i^Ak )  according  to  Equation  (10) 

if  (z(a*)<z(a*)) 

z(A*)<-z(; i‘) 

endif; 

k<^k  +  1 

end; 

Return  best  AP  locations  A  and  overall  objective  value  Z  ( A*  j 

end; 

Figure  3.  Algorithm  DIRECT  for  SRRA+C. 
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is  based  on  the  number  of  APs  n  and  the  number 
of  variables  necessary  to  define  the  location 
of  each  AP  (in  our  case,  a  pair  of  variables  for 
northing  and  easting).  The  vector  \k  represent¬ 
ing  the  locations  of  all  APs  at  iteration  A:  is  a 
single  point  in  this  multidimensional  solution 
space. 

The  initial  solution  A1  is  chosen  as  the  exact 
center  of  the  solution  space.  During  each  itera¬ 
tion,  DIRECT  divides  the  solution  space  into 
smaller  hyper-rectangles,  based  on  the  unex¬ 
plored  territory  in  the  solution  space  and  pre¬ 
viously  calculated  objective  function  values. 
From  any  incumbent  solution,  the  algorithm 
samples  other  locations  in  the  solution  space 
along  each  dimension,  computing  the  corre¬ 
sponding  objective  function  value  for  each.  The 
largest  hyper-rectangles  are  placed  around 
those  sampled  points  with  more  desirable  ob¬ 
jective  values  because  they  have  more  unex¬ 
plored  territory  and  hence  greater  potential  for 
improvement.  During  each  iteration,  the  next 
candidate  solution  \k  depends  on  the  previously 
stored  solution  values  and  the  size  of  their  asso¬ 
ciated  hyper-rectangles.  See  Nicholas  (2009,  pp. 
49-63)  for  more  details  on  our  implementation 
of  DIRECT. 

The  DIRECT  algorithm  is  guaranteed  to 
converge  eventually  to  a  global  optimum  if  the 
objective  function  is  continuous  (Jones  et  al. 
1993).  To  meet  this  requirement,  we  use  bilinear 
interpolation  on  the  discrete  elevation  points 
within  our  map  data  to  smooth  the  terrain 
surface. 

One  drawback  of  the  DIRECT  algorithm  is 
that  the  optimality  gap  at  any  iteration  is  not 
known.  Given  the  Lipschitz  constant,  we  can 
calculate  a  lower  bound,  but  in  our  application 
finding  this  constant  would  be  as  difficult  as 
solving  the  problem  to  optimality.  Another 
drawback  is  that  the  rate  of  convergence  cannot 
be  calculated  exactly.  That  is,  while  DIRECT 
is  guaranteed  to  converge,  we  do  not  know  a 
priori  how  long  this  will  take.  In  general,  this 
is  a  serious  concern,  and  it  could  mean  that 
our  algorithm  runs  too  slowly  to  be  of  practical 
help.  Fortunately,  we  observe  in  practice  that 
it  quickly  yields  solutions  that  are  very  good 
when  compared  to  those  obtained  by  simple 
techniques  like  exhaustive  enumeration  over  a 
finite  grid  of  candidate  solutions.  Though  there 


is  no  guarantee  our  technique  will  always  pro¬ 
vide  good  solutions,  it  is  a  significant  advance 
over  the  kind  of  trial  and  error  currently  used 
by  network  designers  in  the  field. 

As  a  crude  means  of  estimating  an  optimal¬ 
ity  gap  for  a  given  solution,  we  can  calculate 
a  lower  bound  to  the  SRRA+C  objective  value 
by  considering  an  idealized  network  with  zero 
coverage  shortfall  and  no  propagation  loss  be¬ 
tween  APs.  Such  an  idealized  network  does 
not  reflect  the  main  tension  in  our  model:  the 
need  to  spread  APs  to  get  good  coverage  versus 
the  need  to  keep  APs  close  together  to  get  good 
backhaul  capacity.  The  laws  of  physics  and  tech¬ 
nical  limitations  of  the  APs  prevent  us  from 
being  able  to  build  this  idealized  network  in 
practice,  so  this  lower  bound  is  weak,  but  the 
value  provides  a  point  of  reference  against  which 
to  gauge  the  progress  of  the  DIRECT  algorithm. 

SRRA+C  Decision  Support  Tool 

We  implement  the  algorithm  DIRECT  for 
SRRA+C  in  a  decision  support  tool  built  us¬ 
ing  Microsoft  Visual  C+  +  .  The  standalone  pro¬ 
gram  displays  terrain  data,  which  can  be  of 
any  scale  and  any  grid-based  format,  such  as 
Universal  Transverse  Mercator  (UTM)  (Defense 
Mapping  Agency  1989).  The  user  can  input  all 
required  variables,  including  drawing  the  de¬ 
sired  coverage  region(s)  directly  on  a  graphical 
display  of  the  terrain.  The  program  supports 
three  different  modes  of  calculating  path  loss 
(inverse-square,  Hata  COST-231,  and  TIREM). 
After  solving  the  problem,  the  final  AP  loca¬ 
tions,  connections  between  APs,  and  coverage 
areas  are  depicted  atop  the  map.  Note  our  im¬ 
plementation  does  not  require  any  third-party 
solvers  or  additional  software. 


ANALYSIS  AND  RESULTS 

As  case  studies  for  our  analysis,  we  con¬ 
sider  two  different  operating  regions  in  Fort 
Ord,  California.  These  regions  are  respectively 
45  acres  and  145  acres,  and  the  terrain  consists 
of  gently  rolling  hills,  grass  fields,  pavement, 
and  some  buildings.  We  use  terrain  information 
from  the  United  States  Geological  Survey  (USGS) 
via  MapMart  (http:/ / www.mapmart.com). 
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Case  Study  1:  3-node  network 
on  45  acres  in  Fort  Ord,  California 


We  begin  with  a  45-acre  region  (see  Figure  4), 
which  we  divide  into  a  grid  of  size  j R\  =  65  X  33  = 
2,145  rectangles.  We  place  the  HQ  AP  directly 
in  the  middle  of  the  operating  region  and  solve 
SRRA+C  for  a  network  of  three  APs  using 
complete  enumeration.  Enumeration  guarantees 
that  we  will  find  the  optimal  solution  to  this 
restriction  of  the  original  problem  (in  discrete 
space),  but  the  method  is  extremely  resource¬ 
intensive  as  the  number  of  unique  solutions 
grows  exponentially  in  both  locations  |R|  and 
number  of  APs  n.  Recall  the  HQ  node  location 
is  determined  in  advance,  so  building  a  network 
of  n  nodes  means  that  we  need  to  locate  n  —  1 


APs.  We  enumerate  all 


2,299,400 


unique  solutions  in  eight  hours,  38  minutes, 
and  55  seconds.  Figure  5  ranks  the  solutions  by 
their  objective  value. 

Figure  5  shows  that  a  small  number  of  solu¬ 
tions  are  significantly  better  than  the  others. 
Specifically,  about  0.29%  are  within  50%  of 
the  best  solution  obtained,  1,834.8.  Again,  this 
solution  value  has  no  absolute  physical  inter¬ 
pretation  and  serves  only  to  allow  relative  com¬ 
parison.  However,  this  shows  that  the  best 
solutions  are  rare  among  the  possible  solutions. 

We  then  run  DIRECT  on  this  problem  (see 
Figure  6)  and  stop  it  at  each  iteration  1,  2,  ...,  10 


to  record  the  current  overall  objective  value. 
These  values  are  compared  to  our  lower  bound 
for  this  network  (dashed  line)  and  the  best 
value  found  by  enumeration  (solid  line). 

In  only  2.05  seconds,  the  DIRECT  algo¬ 
rithm  obtains  a  better  solution  than  was  found 
by  enumeration  in  8  hours  and  39  minutes. 
This  is  possible  because  enumeration  considers 
only  discrete  locations  within  the  operating  re¬ 
gion,  while  DIRECT  places  APs  continuously 
at  any  location  in  the  region.  Because  it  is  con¬ 
tinuous,  DIRECT  is  guaranteed  to  find  a  solu¬ 
tion  at  least  as  good  as  enumeration  as  the 
number  of  iterations  approaches  infinity.  In 
this  limit,  DIRECT  will  almost  certainly  find  a 
better  one,  because  it  considers  an  increasingly 
dense  subset  of  solutions  and  will  eventually 
sample  to  within  an  arbitrary  distance  of  any 
point  (to  include  the  global  optimum)  within 
the  solution  space. 

In  practice,  networks  with  more  APs  have 
more  dimensions  and  a  corresponding  greater 
number  of  sub-hyper-rectangles.  Hence  for  larger 
networks  the  DIRECT  algorithm  must  be  run 
for  more  iterations  to  find  good  solutions. 
There  is  no  clear  way  to  avoid  this  when  com¬ 
paring  networks  of  varying  numbers  of  APs. 
However,  by  running  it  until  we  reach  the  com¬ 
putational  limit  of  our  implementation  of  the 
algorithm,  we  provide  the  best  answer  possi¬ 
ble.  In  this  case,  running  DIRECT  for  up  to  15 
iterations  does  not  result  in  a  better  solution. 


Elevation  Profile  of  Fort  Ord  Test  Site 


X  Position  (Easting) 


Elevation  (meters) 


Figure  4.  Elevation  profile  of  45-acre  Ford  Ord  region. 
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Figure  5.  Rank-ordered  solutions  for  a  3-node  network  in  45-acre  Fort  Ord  region  obtained  using  enumeration. 


Case  Study  2:  5-node  network 
on  145  acres  in  Fort  Ord,  CA 

We  next  examine  a  network  of  five  APs 
over  a  145-acre  section  of  Fort  Ord.  This  net¬ 
work  is  too  large  to  consider  using  enumera¬ 
tion,  so  we  solve  only  using  DIRECT.  We  again 
stop  the  algorithm  at  each  iteration  1,  2,  ...,  30 
to  record  the  current  overall  objective  value. 
Figure  7  presents  the  results. 


The  best  solution  value  is  computed  by 
DIRECT  after  30  iterations  (796.02)  and  is  within 
an  order  of  magnitude  of  the  lower  bound 
(71.413).  Although  we  cannot  certify  optimality, 
this  boundary  places  a  limit  on  how  far  we  are 
from  the  global  optimum.  Here,  we  are  willing 
to  exchange  optimality  for  speed:  we  obtain  this 
solution  after  22  iterations  of  DIRECT  in  slightly 
more  than  17  seconds;  to  solve  this  problem 
using  our  enumeration  technique  would  take 


Comparison  of  DIRECT  and  Enumeration  Solutions  to  Theoretical  Lower  Bound 


Figure  6.  Comparison  of  solutions  obtained  by  enumeration  and  DIRECT  on  actual  terrain. 
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Comparison  of  DIRECT  Solutions  to  Theoretical  Lower  Bound 
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Figure  7.  Example  comparing  DIRECT  solutions  to  a  lower  bound. 


more  than  104  years.  Of  course,  making  this 
tradeoff  makes  sense  only  if  the  solutions  we 
obtain  are  useful. 

SRRA  +  C  Validation:  Does  It  Work 
in  Practice? 

The  validity  of  solutions  to  the  SRRA+C 
formulation,  specifically  our  evaluation  of  net¬ 
work  performance,  rests  on  the  accuracy  of 
radio  wave  propagation  and  network  through¬ 
put  predictions.  To  validate  our  use  of  TIREM 
to  predict  received  signal  strength  over  real  ter¬ 
rain,  we  conducted  several  field  experiments 
in  collaboration  with  the  Hastily-Formed  Net¬ 
works  (HFN)  Research  Group  at  the  Naval  Post¬ 
graduate  School. 

First,  we  conducted  two  point-to-point 
network  tests  in  Ford  Ord  over  generally  flat 
terrain  consisting  of  pavement  and  packed 
gravel,  with  no  trees  or  other  obstructions  to 
the  line  of  sight  (LOS)  path.  We  deployed 
Cisco  APlOOO-series  Aironet  WMN  access  points 
(Cisco  2009)  broadcasting  at  the  5.8  GHz  and 
2.4  GHz  operating  frequencies  (Figure  8  left  and 
right,  respectively)  and  used  a  laptop  computer 
with  an  internal  Intel  wireless  transceiver  to 
measure  the  received  signal  strength.  We  com¬ 
pared  these  observations  with  TIREM  predic¬ 
tions  at  ranges  from  0  to  465  meters. 

We  observe  that  the  measured  values  are 
reasonably  close  to  the  predicted  ones,  with  the 
exception  of  measurements  at  2.4  GHz  from 
300  meters  and  beyond.  This  is  likely  due  to  the 


presence  of  nearby  buildings  at  that  end  of  our 
testing  range:  the  radio  waves  may  have  re¬ 
flected  off  the  buildings  and  provided  a  stronger 
signal  than  would  have  otherwise  been  re¬ 
ceived.  This  suggests  a  potential  weakness  in 
the  use  of  TIREM  for  urban  environments.  Over¬ 
all,  however,  these  two  tests  demonstrate  that 
TIREM  is  capable  of  making  very  reasonable 
received  signal  strength  predictions  using  our 
testing  equipment  in  a  real-world  environment. 
See  Nicholas  (2009,  pp.  69-75)  for  detailed  testing 
results.  In  particular,  these  results  are  relatively 
accurate  over  ranges  (i.e.,  up  to  250  m)  consistent 
with  having  five  APs  located  uniformly  across 
a  145-acre  area. 

However,  the  most  important  form  of  de¬ 
sign  validation  is  whether  you  can  actually  build 
and  operate  the  system  of  interest.  To  determine 
if  the  topologies  obtained  from  SRRA+C  actu¬ 
ally  function  in  realistic  environments,  we  con¬ 
ducted  another  network  field  test  in  the  45-acre 
rectangular  region  in  Fort  Ord,  mentioned  ear¬ 
lier.  Elevation  here  ranges  from  143  to  226  feet 
above  sea  level.  The  left  side  of  Figure  9  is  a 
Google  Maps  image  of  the  operating  region. 

We  used  the  SRRA+C  Decision  Support 
Tool  to  generate  network  topologies  of  three, 
four,  and  five  APs.  We  deployed  the  same  Cisco 
APlOOO-series  Aironet  WMN  APs,  which  use 
the  802.11b/g  protocol  (2.4  GHz)  to  provide  cli¬ 
ent  coverage  and  the  802.11a  protocol  (5.8  GHz) 
to  provide  the  backhaul  network.  We  posi¬ 
tioned  each  AP  using  a  Global  Positioning  Sys¬ 
tem  (GPS)  device,  and  placed  atop  a  two  meter 
mast.  The  right  side  of  Figure  9  depicts  our  setup. 
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Figure  8.  Validation  of  TIREM  as  a  model  for  received  signal  strength.  At  5.8  GHz  (left)  and  2.4  GHz  (right),  we 
observe  that  our  implementation  of  TIREM  reasonably  predicts  actual  received  signal  strength  in  point-to-point 
field  tests. 


When  required,  we  made  small  adjustments 
(within  five  meters)  to  AP  position  to  avoid 
transmitting  directly  into  trees,  etc.  We  used 
small  portable  generators  and  car  batteries  as 
power  sources. 

Figure  10  displays  the  topologies  prescribed 
by  the  SRRA+C  Decision  Support  Tool  for  the 
4-node  and  5-node  networks.  After  building 
each  network,  we  used  a  laptop  to  measure  the 
actual  throughput  from  each  AP  to  the  HQ 
node  ( d ).  The  results  appear  alongside  each 
node  in  Figure  10. 

The  network  topologies  in  Figure  10  worked 
similarly  to  our  theoretical  predictions.  The 
backhaul  network  traffic  from  each  AP  to  the 
HQ  node  followed  exactly  the  routes  predicted 
by  SRRA+C.  The  actual  throughput  from  each 


AP  to  the  HQ  node  was  within  an  order  of 
magnitude  of  our  predictions,  but  typically 
exceeded  them.  While  the  throughput  values 
predicted  by  SRRA+C  were  based  on  the  as¬ 
sumption  that  all  APs  were  transmitting  con¬ 
currently,  in  our  experiments  we  were  limited 
to  measuring  the  throughput  from  each  AP  in¬ 
dividually.  Because  a  lone  transmitting  node 
does  not  have  to  compete  for  backhaul  network 
resources,  it  is  natural  that  observed  through¬ 
put  should  exceed  prediction.  Given  our  sim¬ 
plifying  assumptions  about  the  way  the  APs 
work  (e.g.,  ignoring  the  signal  processing  tech¬ 
nologies  that  try  to  compensate  for  environ¬ 
mental  conditions),  the  level  of  correspondence 
between  the  predicted  and  observed  values 
exceeded  our  expectations. 


Figure  9.  Field  testing  in  Fort  Ord,  California.  Left:  Aerial  view  of  45-acre  Fort  Ord  operating  region  [image 
courtesy  Google  Maps  (http:/ / maps.google.com)].  Right:  Wireless  mesh  access  point,  mast,  and  portable  generator 
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Figure  10.  Network  topologies  and  their  measured  performance  during  field  tests.  Circles  denote  position  of 
each  AP  during  field  tests  of  4-node  (left)  and  5-node  (right)  networks.  For  each  node,  we  report  the  actual 
(and  predicted)  throughput  from  that  node  to  the  HQ  node,  denoted  d.  Arrows  represent  the  direction  of  actual 
point-to-point  wireless  transmissions  along  the  backhaul  network,  which  exactly  matched  the  predicted  flows. 
For  the  5-node  network,  SRRA+C  predicted  that  node  5  should  split  its  transmission  in  two  directions;  we  ob¬ 
served  that  this  node  alternatively  sent  flow  in  two  directions. 


Not  all  of  our  field  tests  were  successful. 
In  the  three-AP  scenario  (not  shown),  a  node 
was  unable  to  connect  with  the  HQ  node  and 
the  throughputs  were  lower  than  expected.  We 
attribute  this  to  the  presence  of  intervening  fo¬ 
liage  and  man-made  obstructions,  which  are  not 
considered  by  the  TIREM  model.  Nonetheless, 
our  field  testing  results  show  that  SRRA+C 
solved  with  DIRECT  can  quickly  provide  work¬ 
ing  network  designs  with  no  guesswork.  See 
Nicholas  (2009,  pp.  92-102)  for  a  detailed  de¬ 
scription  of  these  tests  and  results. 

How  Many  APs  Do  We  Need? 

The  SRRA+C  problem  allows  us  to  deter¬ 
mine  where  to  place  a  known  number  of  APs 
within  a  given  operating  region.  But  an  equally 
important  question  is:  How  many  APs  do  we 
need ? 

To  address  this  question,  we  again  consider 
the  45-acre  section  of  Fort  Ord  and  we  use 
DIRECT  to  solve  for  networks  of  2,  3,  ...,  8 
APs  using  weight  w  =  1,  running  the  algorithm 
until  we  reach  the  computational  limit  of  32 
hyper-rectangle  divisions  (the  maximum  possi¬ 
ble  using  double-precision  floating  point  num¬ 
bers).  We  plot  the  best  solution  obtained  for 
each  network  as  a  function  of  the  two  competing 
terms  in  the  overall  objective  function  (Equation 
10).  The  information  presented  in  such  a  plot 
serves  two  purposes.  First,  it  provides  informa¬ 


tion  on  the  relative  "goodness"  of  a  particular 
solution  by  comparing  it  to  the  bounds  of  net¬ 
work  flow  and  coverage  shortfall.  It  also 
provides  a  quantification  of  the  value  of  addi¬ 
tional  APs  and  the  appropriate  number  of  APs 
for  a  particular  scenario. 

Figure  11  presents  the  best  solutions  to  net¬ 
works  having  two  to  eight  APs  in  the  45-acre 
Fort  Ord  region  and  illustrates  the  relative 
contribution  of  client  coverage  and  network 
flow  to  the  overall  objective  value.  We  observe 
a  decreasing  gain  in  coverage  with  more  than 
four  APs,  which  we  interpret  to  mean  that  four 
APs  are  sufficient  to  provide  coverage  to  the 
majority  of  that  region.  As  expected,  additional 
APs  provide  greater  network  flow. 

Observe  the  main  improvement  from  the 
optimal  2-node  network  to  the  optimal  3-node 
network  is  a  decrease  in  coverage  shortfall 
(and  in  fact  there  is  a  slight  decrease  in  delivered 
network  throughput).  The  optimal  4-node  net¬ 
work  provides  significantly  higher  delivered 
network  throughput  while  also  reducing  the 
coverage  shortfall.  For  networks  of  size  five 
and  larger,  additional  APs  primarily  increase 
delivered  network  throughput.  Interestingly,  it 
appears  that  only  after  achieving  adequate  cli¬ 
ent  coverage  does  it  make  sense  to  use  addi¬ 
tional  APs  to  enable  greater  network  flow. 

This  chart  helps  answer  the  questions 
"How  good  is  our  network?"  and  "How  much 
better  can  we  do?"  This  information  is  of 


Military  Operations  Research,  V17  N4  2012 


Page  81 


FAST,  EFFECTIVE  TRANSMITTER  PLACEMENT  IN  WIRELESS  MESH  NETWORKS 


0  1000  2000  3000  Less  coverage 

Total  coverage 

Coverage  Shortfall 

Figure  11.  Example  bi-objective  plot  of  SRRA+C  for  the  45-acre  Fort  Ord  region.  Each  point  represents  our  best 
solution  to  Z(A)  for  the  specified  number  of  APs,  solved  with  equal  objective  weight  given  to  the  coverage  and 
flow  terms  (w  =  1).  Beyond  four  APs,  we  observe  decreasing  benefit  in  coverage  shortfall  and  increasing  benefit 
in  delivered  network  flow. 


considerable  value  to  a  decision-maker,  and  the 
speed  of  our  SRRA+C  algorithm  makes  this 
analysis  possible. 


CONCLUSION  AND 
RECOMMENDATIONS 
FOR  FUTURE  WORK 

In  this  paper,  we  present  a  technique  for  de¬ 
signing  WMNs  that  maximizes  client  coverage 
area  while  considering  network  service  and  ra¬ 
dio  propagation  over  terrain.  Our  SRRA+C  for¬ 
mulation  does  not  require  information  about 
device-specific  characteristics  such  as  radio 
wave  modulation  scheme  or  network  routing 
protocol.  Although  these  simplifications  may 
reduce  the  predictive  power  of  the  formula¬ 
tion,  it  is  possible  to  add  such  considerations. 
Further,  this  general  approach  makes  it  very 
easy  to  quickly  model  networks  of  diverse  de¬ 
vices  and  capabilities. 

Using  our  decision  support  tool,  we  com¬ 
pare  the  performance  of  enumeration  and  the 
DIRECT  algorithm  in  solving  the  SRRA+C  prob¬ 


lem.  We  show  that  the  DIRECT  algorithm  can 
find  good  solutions  much  faster  than  complete 
enumeration,  and  is  capable  of  finding  better  so¬ 
lutions  than  enumeration  because  DIRECT  con¬ 
siders  a  continuous  solution  space.  DIRECT  is 
guaranteed  to  find  the  optimal  solution  to  the 
SRRA+C  problem  as  the  number  of  itera¬ 
tions  goes  to  infinity,  and  each  iteration  of 
DIRECT  provides  a  solution  at  least  as  good 
as  the  previous. 

Although  we  cannot  provide  certificates 
that  guarantee  the  optimality  of  any  SRRA+C 
solution,  proof  of  optimality  in  this  domain  is 
not  necessary  and  in  fact  wastes  critical  time. 
In  practice,  operators  need  good  working  WMN 
topologies  quickly,  and  our  technique  provides 
exactly  this. 

Our  techniques  and  associated  decision  sup¬ 
port  tool  can  be  used  by  HA/ DR  personnel  and 
combat  communications  planners  to  quickly  de¬ 
sign  WMNs  to  support  their  specific  operations. 
The  software  runs  on  a  laptop  computer,  does 
not  require  any  additional  software  or  solver 
licenses,  accepts  map  data  in  a  generic  file  for¬ 
mat  that  is  widely  available  on  the  Internet, 
and  creates  network  topologies  for  virtually 
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any  type  of  terrain  and  mesh  AP  device.  This 
technique  requires  very  little  technical  expertise 
and  no  guesswork. 

Future  work  could  investigate  the  use  of 
higher-fidelity  radio  propagation  and  network 
traffic  models,  and  compare  the  benefit  of 
increased  accuracy  to  any  additional  computa¬ 
tional  workload.  In  particular,  we  anticipate  the 
need  to  investigate  radio  propagation  models 
that  consider  vegetation  and  man-made  obstruc¬ 
tions,  which  are  critical  in  calculating  radio  prop¬ 
agation  in  urban  environments. 

We  have  shown  that  complete  enumeration 
has  limited  usefulness  as  a  solution  technique 
for  comparison  to  DIRECT.  Future  work  could 
compare  DIRECT  to  other  sampling  algorithms 
or  heuristic  approaches.  Future  work  could 
also  explore  the  use  of  DIRECT  in  a  parallel- 
or  multiple-processor  architecture  as  DIRECT 
naturally  lends  itself  to  this  type  of  computa¬ 
tion  (see,  e.g..  He  et  al.  2004). 

We  are  currently  exploring  the  use  of 
SRRA+C  with  DIRECT  to  build  larger  net¬ 
works  (of  100  APs  and  more)  and  networks 
with  more  than  one  traffic  destination.  Assessing 
the  impact  of  disruptions  (accidental  or  inten¬ 
tional)  on  the  design  of  WMNs  (e.g.,  Grotschel 
et  al.  1995,  and  Shankar  2008)  is  a  topic  that  also 
warrants  additional  study. 
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